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Abstract 

Moncher is a Monte Carlo event generator for simulation of single and double 
charge exchange reactions in proton-proton collisions at energies from 0.9 to 14 TeV. 
Such reactions, pp — > n + X and pp — > n + X + n, are characterized by leading 
neutron production. They are dominated by 7r + exchange and could provide us with 
more information about total and elastic ir + p and 7i + tt + cross sections and parton 
distributions in pions in the still unexplored kinematical region. 
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1 Introduction 

In the paper we present a new Monte-Carlo event generator Moncher. The generator 
is devoted to the simulation of single and double charge exchange reactions in proton- 
proton collisions at energies from 0.9 to 14 TeV. This region of energies covers the present 
capabilities of the LHC Charge exchange reactions, pp — > n + X and pp — > n + X + n, are 
characterized by the leading neutron production. They can be studied with LHC detec- 
tors incorporated with forward neutron calorimeters like the ZDC (Zero Degree Calorime- 
ter) [I] in the CMS 0. 

Reactions with the leading neutron production are dominated by 7r + exchange [3]- [6]. 
At the LHC they could provide us with information about ir + p and tt + tt + interactions in 
the region of energies 1-5 TeV in the c.m.s. Using indirect methods [7]- [8] we could extract 
total and elastic 7t + p and tt + tt + cross sections at these energies. It is worth mentioning 
that the total cross-section of 7r + p interaction is measured only at energies up to 25 GeV 
by direct methods in the fixed target experiments [9] and total and elastic cross sections of 
7r + 7r + interactions are extracted from the data at energies 1.5-18.4 GeV only (see Ref. |10j - 
[12]). Moreover, a study of charge exchange reactions with hard scattering 7c + p and tt + tt + 
followed by dijet production at the LHC, could provide us with parton distributions in the 
pion in the unexplored kinematical domain. So, we had weighty motivations to develop a 
model and to create a generator for charge exchange simulation which could be used at 
high energies of the LHC. 

An important point is that at high energies we have to take into account effects of 
soft rescattering which can be calculated as corrections to the Born approximation. In 
the calculations of such absorptive effects we use the Regge-eikonal approach [T3]. For 
tt + p and 7r + 7r + interactions several models which predict different cross sections have 
been applied. In addition to the dominant 7r + exchange we have calculated contributions 
of two other important Reggeons, p + and a^, to the charge exchange cross section p3] 
and implemented both Reggeons to the generation. Pythia 6.4 (15] is used as a basic 
generator for Moncher. Moncher has the same format of events, parameters and 
common blocks as Pythia . PYTHlAsubroutines are used also for the simulation of 7i + p 
and 7T + 7r + interactions and for the subsequent hadronization and decays. 



2 Physics Overview 
2.1 Single Pion Exchange 

The diagram of the Single pion Exchange (S7rE) process p + p — > n + X is presented in 
Fig. [T^l. The momenta are pi, P2, p n , Px respectively. In the center-of-mass frame these 
can be represented as follows (boldface letters denote transverse momenta): 

With this notation, the momentum of the 7T is 



and 



Pn=Pl-p7T, (3) 

Px = M\ (4) 
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Figure 1: Amplitudes of the processes: a) p + p — > n + X (S7rE), b) p + p — > n + X + n 
(D-7rE). 5" and 5*2 represent soft rescattering corrections. 
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As a Born approximation for tt exchange we use the familiar triple-Regge formula. 
This formula can be rewritten as 



da 



X,SirE 
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d^dtd^x 16vr 2 (t - m 2 ) 

dax^+pi^s) 
x ^ S(s/s ,Z,t), 
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where $x is the phase space for the system X produced in the 7r + p scattering, the pion 
trajectory is a w (t) = a' n (t — m 2 ). The slope a' ~ 0.9 GeV -2 , £ = 1 — x^, were xl is 
the fraction of the initial proton longitudinal momentum carried by the neutron, and 
G 2 0pp /(47r) = G 2 +pn /(Sii) = 13.75 [TBI [17] . The form factor F (t) is usually expressed as 
an exponential 

F (t) = exp(fa), (9) 

where, from recent data [18], [19], we expect b ~ 0.3 GeV -2 . We are interested in the 
kinematical range 

0.01 GeV 2 < \t\ < 0.5 GeV 2 , £ < 0.4, (10) 

where formula (JHJ) dominates according to [20j and [21]. At high energies we can use any 
adequate parametrizations of different n + p cross-sections. 



2.1.1 Absorptive corrections 

The suppression factor S arises from absorptive corrections [3]. We estimate absorption 
in the initial state for inclusive reactions and for both initial and final states in exclusive 
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exchanges. For this task we use our model with 3 Pomeron trajectories 

a IPl (t) - 1 = (0.0578 ± 0.002) + (0.5596 ± 0.0078)t , 

a IP2 (t) - 1 = (0.1669 ± 0.0012) + (0.2733 ± 0.0056)t , (11) 

a IP3 (t) - 1 = (0.2032 ± 0.0041) + (0.0937 ± 0.0029)t . 

These trajectories are the result of a 20 parameter fit of the total and differential cross- 
sections in the region 

0.01 GeV 2 < |t| < 14 GeV 2 , 8 GeV < < 1800 GeV. 

Although the \ 2 jd.o.j. = 2.74 is rather large, the model gives good predictions for the 
elastic scattering (especially in the low-t region with x 2 /d.o.f. ~ 1). 

We use the procedure described in [5], [6] to estimate the absorptive corrections. With 
an effective factorized form of (see hereunder) expression fll2p used for convenience, we 
obtain: 

== He + ^s^^^^a^ s), (13) 

_ mie\Ms/s , g, q 2 )| 2 + q 2 |<fr a (*/*o, g, q 2 )| 2 

(m 2 £ 2 + q 2 )|<M£,q 2 )| 2 ' 1 ' 

The functions $o and $ s arise from different spin contributions to the amplitude 

A p ^ n = (m P i d 3 • $ + q & ■ $ s ) (15) 

a/i — 4 

and both are equal to in the Born approximation. Here &\ are Pauli matrices and 
ty n , \l/p are neutron and proton spinors. All the above functions can be calculated by the 
following set of formulae: 

-■"<« 1 1 q-0, (16) 
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Ci 


53.0 ±0.8 


9.68 ±0.16 


1.67 ±0.07 


rf (GeV- 2 ) 


6.3096 ±0.2522 


3.1097 ±0.1817 


2.4771 ± 0.0964 



Table 1: Parameters of the model. 



V(b)=exp(-n el (s/s ,b)) 

3 

i=l 
2d 



a- 



B, 



167T Bi \Sq 



— e 2 ) exp 



45, 



a' m ln(^f) + ^. 



(23) 
(24) 

(25) 
(26) 



The values of parameters Cj and rf are derived in (TlTI) and listed in Table [TJ Figs. [2] 
demonstrate function S(s/s , £, %) calculated for two values of energies a) y/s = 62.7 GeV 
and b) y/s = 10 TeV for different £ values: £ = 0.3 (dotted), £ = 0.1 (dashed) and £ = 10~ 4 
(solid). 
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Figure 2: Function S(s/s ,£,,q t ) at a) y/s = 62.7 GeV and b) y/s 
£ values: f = 0.3 (dotted), f = 0.1 (dashed) and f = 10" 4 (solid). 
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2.1.2 Parametrization of ir + p cross section 

In the present version of generator we use 4 parametrizations for 7r + p cross section. 

The Donnachie-Landshoff (DL) parametrization |22j: 

a%{s) = 13.63 s - 0808 + 25.56 s" ' 4525 , (mb). (27) 
The COMPETE parametrization [23]: 

a%{s) = Z„ p + Bin 2 (-) + (Y + s a + - Y_s a ~) /s, (mb). (28) 
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Table 2: Parameters of the model 12 



Z wp = 21.23 ± 0.33 mb, 

B = 0.3152 ±0.0095 mb, 

s = 34 ±5.4 GeV 2 , 

Y + = 17.8 ± 1.10, a + = 0.533 ± 0.015, 

F_ = 5.72 ± 0.16, a_ = 0.4602 ± 0.0064. 

In the next two parametrizations total cross-section can be obtained thorough the 
optical theorem 
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The Bourrely-Soffer-Wu (BSW) parametrization [24J : 
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For the ir + p we have i = p in ( 131]) and 
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B p = b p + a' (0) In — , b p = 4.2704, 

So 

a p (t) = 0.3202 + t, C p = 4.1624, 
where values of parameters are listed in Table [2j 
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The Godizov-Petrov (GP) parametrization |25j,[26j. 

In this parametrization the scattering amplitude is represented in the usual eikonal form 

e 2iS(s,b) _ i 



T(s,b) 



2i 



(37) 



(here T(s, b) is the amplitude in the impact parameter b space, s is the invariant mass 
squared of colliding particles and S(s, b) is the eikonal function). Amplitudes in the impact 
parameter space and momentum one are related thorough the Fourier-Bessel transforms 
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Eikonal function in the momentum space is 

5{s,t) = S P (s,t) + S f (s,t) = 

I ,rMt)-l)\ /.N-M 



Z + tE 



Vs 



+ 



/. 7r(a f (t) W') 

The parametrization for the pomeron residue is 

/3 p (t) = S P e^*(l + di t + d 2 1 2 + d 3 1 3 + di t 4 ) , 

which is approximately (at low values of d\, d 2 , d 3 <i 4 ) an exponential at low t values. 
Residues of secondary reggeons we set as exponentials: 

j 8 / (t) = S / e 6 /*. (42) 
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b P 


2.4 GeV" 2 


bf 


4.7 GeV" 2 


K 


5.6 GeV- 2 
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a' P (0) 


0.28 GeV" 2 


«f(0) 


0.63 GeV" 2 




0.07 GeV" 2 



Table 3: Values of parameters of the model [25], [26] for pp scattering. 
Phenomenological parametrization for the "soft" pomeron trajectory is set to 
a P (t) = 1 +pi 
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l-p 2 t(&TCtg(p 3 -p 2 t) _ } 

Trajectories of secondary reggeons f 2 and u are parametrized by functions 
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is the one-loop analytic QCD running coupling [27], rif = 3 is the number of flavours, 
A = A (3) = 0.346 GeV [28]. Parameters Cf i c u > are rather small to spoil the asymptotic 
behaviour of secondary trajectories in the perturbative domain. Residues for tttc, Tip and 
pp are assumed to be 



PF{t) 



0?(t)0?{t) 



m p (t) 



(47) 



Parameters of the model are listed in Tables [3j 
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2.2 Double Pion Exchange 

The diagram of the Double Pion Exchange (D7rE) process p + p — > n + X + n is presented 
in Fig. [Tb. The momenta are pi, p 2 , p ni , Px, Pn 2 respectively. In the center-of-mass frame 
these can be represented as follows 



Pni Pi Pni , 

Px — M 2 = 66^^ 



+ t 2 + 2(m 2 - mi)) ■ (£ 2 (6 + 6) + 



t 1 +t 2 + 2(m 2 -m 2 ) 
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66«, (50) 
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The cross-section can be evaluated as follows: 
da = S 2 (s/s ,£ 1>2 ,<ll 2 )d<TQ, 
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For low tj the function is approximately equal to 

F(6,6) = S2(sAo,£i,£ 2 ,o,o)~ 
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^(5/50,6,0)+ y/^s/so, 6,0) - yJS(s/s ,^,0)S(s/s Q ,^,0)) . (5 



Figs. |3] demonstrates 2D projections of function S 2 (s/so, 6,2, | 2 1 ) an d function ^(6,6) 
at ^ = 10 TeV. 

To obtain cross-sections we use parametrizations described in the subsec- 

tion 12.1.21 with the following approximations: 
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Figure 3: The function S 2 (s/s , £i j2 , |q x 2 |) at y/s = 10 TeV: a) for fixed £i )2 = 0.01 b) for 
fixed |q x 2 | ~ 0. c) The function F(£i,£ 2 ) at a/s = 10 TeV. 
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for the BSW one, which should be substituted to ( 132|) . and equations ( H6]) .( 14"T|) for the 
GP one. 

2.3 Relative contributions of 7r, p and a2 reggeons. 



For p and a 2 contributions formulae are similar to ones described in the chapters (12. ip 
and (ED- 
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Here % = 8 is the ratio of spin-flip to nonflip amplitude, aja(t) — 0.5 + 0.9t and param- 
eters for p, a>2 mesons are [29] 



V P = -i+l, ^a 2 =i + l, 

6 P = 2 GeV~ 2 , 6 a2 = 1 GeV~ 2 , 



8tt 



8tt 



(67) 
(68) 

(69) 



Rescattering corrections Sr and Sr^ are calculated by the method used in [7], [8]. Basic 
assumptions in our calculations are: 



p p, p a 2 and a 2 a 2 contributions are small; 
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• interference terms of the type T^Ts^e, T^^T^^i^ are small [19], 11,11' = 
ii, p, a 2 , 1R 7^ IR', where T are amplitudes of the corresponding processes; 

• approximate relations %+ p ~ a n + p , (7^+^+ ~ a 7r + 7r + [T9] . 

Figs. H] demonstrates 3D plots for cross sections of the Single and Double Reggeon 
Exchange reactions for the different reggeons (S7rE a), SpE+Sa 2 E b), D7rE c) and 
Dp7rE+Da 2 vrE d)) at yfs = 7 TeV. 




Figure 4: Cross-sections in mb ■ cm 1 at -Js = 7 TeV for: a) S7rE; b) SpE+Sa 2 E; c) 
D7rE; d) Dp7rE+Da 2 7rE. r is the transverse distance from the beam. 
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3 Program Overview 



The kinematics of SHE and DIRE reactions, 

pp^n + (TRp) -> n + X (70) 

and 

pp -)■ ra + (RE) + n-»n + X + n, (71) 

are defined by the relative energy loss of £ n and the square of the transverse momentum 
t n of the leading neutron. The vertex pEn is generated on the basis of the models 
described above. The differential cross sections for the generated neutron and reggeon 
are calculated according the selected models for absorptive corrections and for JRp (1R1R) 
interactions. Then, PYTHIA 6.420 [15] is called for the Hp — > X generation in the case 
of SIRE and 1R1R — >■ X generation in the case of DIRE. Parameters of the all generated 
particles, including beam protons, leading neutrons, reggeons and X, products of IRp 
(1R1R) interaction, are stored in PYTHlAcommon blocks. 

3.1 Main Subroutines 



SUBROUTINE MONINIT 



Purpose: to initialize the generation procedure. In particulary, 

- to show program title; 

- to read control parameters from the file moncher.par; 

- to set default Moncher parameters; 

- to initilize PYTHIA; 

- to initilize LHE format output. 

Status of call : should be called obligatory, one time, in the begining of the main 

program before calling of MONEVEN. 
Calling by: main program 

Calling of : MONTITL, MONPARA, MONMBDF, MONUPIN, PYINIT 



SUBROUTINE MONTITL 



Purpose: to print title of Moncher on the screen. Namely, 



* * 

* MON-te-carlo generator for * 

* CH-arge " * 

* E-xchange * 

* R-eactions * 
^ ^ 

* Version 1.1.0. (12/03/2011) * 

* * 

* * 

* * 

* R.Ryutin.A.Sobol.V.Petrov * 

* (IHEP,Protvino) * 

* * 



Status of call : should be called OBLIGATORY. 
Calling by: MONINIT 
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SUBROUTINE MONPARA 



Purpose: to read control parameters for the Moncher and Pythia generation from 
the file moncher. par Control parameters for MONCHER are called MONPAR, 
they are stored to the common block /MONGLPA/. Any PYTHIAparameters 
can be defined for PYTHlAcommon blocks 

/PYJETS/ , /PYDAT1/ , /PYDAT2/ , /PYDAT3/ , /PYDAT4/ , /PYDATR/ , /PYSUBS/ , 
/PYPARS/ , /PYINT1/ , /PYINT2/ , /PYINT3/ , /PYINT4/ , /PYINT5/ , /PYINT6/ , 
/PYINT7/,/PYINT8/,/PYMSSM/,/PYMSRV/,/PYTCSM/,/PYPUED/ (see [15]). 

Status of call : can be called if you like to define some control parameters from the 
moncher .par. By default, MONCHER initilizes a generation of minimum bias 
events by Pythia with some default parameters. 

Calling by: MONINIT 

Calling of: MONGIVE 



SUBROUTINE MONMBDF 



Purpose: to define default parameters for the generation. By default, Moncher and 
Pythia parameters are defined to generate 10 minimum bias events at c.m.s. 
energy 7 TeV. 

Status of call : to be called at initialization. Default parameters are redefined by the 

call of the MONPARA reading parameters from the file moncher . par. 
Calling by: MONINIT 
Calling of: MONGIVE 



SUBROUTINE MONEVEN 



Purpose: call subroutines for the single event generation 

M0NPAR(7)= 1 : call MONSPEG for Single Charge Exchange (SCE) generation. 

M0NPAR(8)= 1 : call MONDPEG for Double Charge Exchange (DCE) generation. 

M0NPAR(7)= and M0NPAR(8)= : call PYEVNT for the PYTHIAevent generation. 
Status of call : should be be called in the user main program, in the cycle of events. 
Calling by: main program 
Calling of : MONSPEG, MONDPEG, PYEVNT 



SUBROUTINE MONSPEG 



Purpose: to generate single SCE event, p b ^ am p b ^ am _>. n(TT^ irt p2 eam ) — > nX, in the follow- 
ing sequence: 

- the vertex p 1 eam n7r+ rf is generated by MONSPEM; 

- Pythia is initialized for the generation of ^tirtP b 2 am interaction; 

- Pythia is called for the generation, hadronization and decays; 

- the PYTHlAoutput is rewriting to include beam protons and neutron 
to the final state of the reaction with the particles from X. 

Simulation of 7r^ irt p b f am interaction is controled by PYTHIAparameters. It can 
be elastic, minimum bias or diffractive interaction. Number of the correspond- 
ing SCE process is equal to the number of the PYTHlAprocess + 500. 

Status of call: called if MONPAR (7) =1. 

Calling by: MONEVEN 

Calling of : MONSPEM, MONSHPY, PYINIT, PY1ENT, PYANGL 
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SUBROUTINE MONDPEG 



Purpose: to generate single DCE event, p b ^ am p h ^ am 



— > n{ix, 



vtrt vtrtJ 



nXn, in the 



following sequence: 

- the vertexes Pi eam "nir^ irt and p^ am rnr^ irt are generated by MONDPEM; 

- Pythia is initialized for the generation of 7r+ rt 7r+ rt interaction; 

- Pythia is called for the generation, hadronization and decays; 

- the PYTHlAoutput is rewriting to include beam protons and neutrons 
to the final state of the reaction with the particles from X. 

Simulation of 7r^ rt 7r^ rt interaction is controled by PYTHlAparameters. It can 
be elastic, minimum bias or diffractive interaction. Number of the correspond- 
ing DCE process is equal to the number of the PYTHlAprocess + 600. 

Status of call: called if M0NPAR(8)=1. 

Calling by: MONEVEN 

Calling of : MONDPEM, MONSHPY, PYINIT, PY1ENT, PYANGL 



SUBROUTINE M0NSPEM(N0,PN,PR,M2) 



Purpose: to generate momentums and energies of neutron n and virtual exchange 
reggeon 1R + in the reaction of Single Charge Exchange: 



p beam p beam _^ n (TR+ p beamj _^ nX 



INTEGER NO 
= 1 
= 2 
= 3 



( input) 

7T + 
P + 

at 



type of exchange reggeon H; 



DOUBLE PRECISION 



Px, 
Pz, 



PN(1) 
PN(2) 
PN(3) 
PN(4) 
PN(5) 
DOUBLE PRECISION 
PR(1) 
PR(2) 
PRO) 
PR(4) 
PR(5) 



DOUBLE PRECISION M2 (output) 
Calling by: MONSPEG 
Calling of: M0NGE2D 



PN(5) (output) : kinematical parameters of the neutron n. 
momentum of neutron in the x direction, in GeV/c. 
momentum of neutron in the y direction, in GeV/c. 
momentum of neutron in the z direction, in GeV/c. 
E, energy of neutron, in GeV. 
m, mass of neutron, in GeV/c 2 . 

PR(5) (output) : kinematical parameters of the reggeon IR + . 
p x , momentum of reggeon in the x direction, in GeV/c. 
p y , momentum of reggeon in the y direction, in GeV/c. 
p z , momentum of reggeon in the z direction, in GeV/c. 
E, energy of reggeon, in GeV. 
m, mass of reggeon, in GeV/c 2 



invariant mass of the system (\R + p beam2 ), in GeV/c . 



SUBROUTINE MONDPEM (NO , PN1 , PN2 , PR1 , PR2 , M2) 



Purpose: to generate momentums and energies of neutrons n and virtual exchange 
reggeons H + in the reaction of Double Charge Exchange: 

p bea mp bea m _^ n (]R+]R+) n n Xn 

INTEGER NO (input) : type of exchange reggeons IR^IR^; 



= 1 

= 2 
= 3 



7T 7T 
7T + p + 



DOUBLE PRECISION PN1 (5) ,PN2 (5) (output) : kinematical parameters of the neu- 
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tron n. 

PN1(1) ,PN2(1) : p x , momentum of neutrons in the x direction, in GeV/c. 
PN1(2) ,PN2(2) : p y , momentum of neutrons in the y direction, in GeV/c. 
PN1(3) ,PN2(3) : p z , momentum of neutrons in the z direction, in GeV/c. 
PN1(4) ,PN2(4) : E, energy of neutrons, in GeV. 
PN1(5) ,PN2(5) : m, mass of neutrons, in GeV/c 2 . 
DOUBLE PRECISION PR1 (5) ,PR2 (5) (output) : kinematical parameters of the 
reggeons H^ 2 - 

PR(1) ,PR2(1) : p x , momentum of reggeons in the x direction, in GeV/c. 

PR(2) ,PR2(2) : p y , momentum of reggeons in the y direction, in GeV/c. 

PR(3) ,PR2(3) : p z , momentum of reggeons in the z direction, in GeV/c. 

PR (4) ,PR2(4) : E, energy of reggeons, in GeV. 

PR(5),PR2(5) : m, mass of reggeons, in GeV/c 2 . 
DOUBLE PRECISION M2 (output) : invariant mass of the system (M|IRJ), in GeV/c 2 . 
Calling by: MONDPEG 
Calling of : M0NGE2D, M0NGE2D4 



SUBROUTINE MONSHPY (NSHIFT) 



Purpose: to shift data of arrays of the PYTHIAcommon block /PYJETS/ for NSHIFT 
positions. It should be done to fill first NSHIFT positions of /PYJETS/ arrays 
by the parameters of the beam protons and neutrons in the final state of 
reaction. 

Calling by: MONSPEG, MONDPEG 



SUBROUTINE MONGIVE(CHIN) 



Purpose: modification of the PYTHIAsubroutine PYGIVE to set the value of any variable 
residing in the commmonblocks PYJETS, PYDAT1, PYDAT2, PYDAT3, PYDAT4, 
PYDATR, PYSUBS, PYPARS, PYINT1, PYINT2, PYINT3, PYINT4, PYINT5, PYINT6, 
PYINT7, PYINT8, PYMSSM, PYMSRV, PYTCSM or MONGLPA. This is done in a more 
controlled fashion than by directly including the common blocks in your pro- 
gram, in that array bounds are checked and the old and new values for the vari- 
able changed are written to the output for reference. In the following example, 
"CALL M0NGIVE('M0NPAR(3) = 14000')", we have changed pp c.m.s. energy to 
14 TeV. More detail explanation see in Ref. (15] for subroutine PYGIVE. 

CHARACTER CHIN*(*) (input) : character expression of length at most 100 characters, 
with requests for variables to be changed. 

Calling by: MONPARA, MONMBDF 



SUBROUTINE MONUPEV 



Purpose: to write information about generated processes to the file moncher . lhe using 
special LHE record format. For more detail information about LHE format 
see Ref. [3D]. 

Status of call: called if M0NPAR(2)=1. 

Calling by: MONINIT 



SUBROUTINE MONUPIN 



Purpose: to save information about all stable particles generated in the event to the file 
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moncher.lhe using special LHE record format. For more detail information 

about LHE format see Ref. [30] . 
Status of call: should be called for each generated event if MONPAR (2) =1. 
Calling by: user main program 

3.2 Auxiliary Subroutines 

These subroutines are used for internal calculations and should not be changed. 



SUBROUTINE M0NGE2D (FF , XI , X2 , Nl , N2 , FF1 , FF2 , FF3 , RG , XG , IG) 



Purpose: to generate two variables according to the 2D distribution from the table. 
DOUBLE PRECISION FF(N1 ,N2) (input) : NlxN2 dimensional interpolation table of 2D 
distribution. 

DOUBLE PRECISION XI (Nl) ,X2(N2)) (.input) : arrays of variables corresponding to the 
table FF. 

INTEGER Nl ,N2 (input) : dimensions of the 2D table. 

DOUBLE PRECISION FFl(Nl) ,FF2(N1) (input) : auxiliary integrated tables for 2D dis- 
tribution. 

DOUBLE PRECISION FF3(2,N1,N2) (input) : auxiliary sums from the table for 2D distri- 
bution. 

DOUBLE PRECISION RG(2) (input) : array for generated random numbers from to 1. 
DOUBLE PRECISION XG(2) (output) : array for generated variables according to the 2D 
distribution. 

INTEGER IG(2) (output) : auxiliary numbers of the nearest to the XG (2) discrete point. 
Calling by: MONSPEM, MONDPEM 



SUBROUTINE M0NG2D4 (FF , XI , X2 , Nl , N2 , FF1 , FF2 , FF3 , I I , XX , RG , XG) 



Purpose: to generate four variables according to the 4D distribution from the table. 

DOUBLE PRECISION FF(N1 , Nl ,N2 ,N2) (input) : NlxNlxN2xN2 dimensional interpola- 
tion table of 4D distribution. 

DOUBLE PRECISION XI (Nl) ,X2(N2)) (input) : arrays of variables corresponding to the 
table FF. 

INTEGER m ,m(input) : dimensions of the 4D table. 

DOUBLE PRECISION FF1 (4,N1 ,N1 ,N2) ,FF2(4,N1 ,N1 ,N2) (input) : auxiliary integrated 

tables for the 4D distribution. 
DOUBLE PRECISION FF3 (8 , Nl ,N1 , N2 ,N2) (input) : auxiliary sums from the table for the 

4D distribution. 

INTEGER 11(2) (input) : auxiliary numbers for multidimensional calculations. 
DOUBLE PRECISION XX(2) (input) : auxiliary points for multidimensional calculations. 
DOUBLE PRECISION RG(4) (input) : array for generated random numbers from to 1. 
DOUBLE PRECISION XG(4) (output) : array for generated variables according to the 4D 

distribution. 
Calling by: MONDPEM 



SUBROUTINE MONCUBI (FF,VS,FUN) 



Purpose: cubic spline interpolation for a function in the variable In s. 
DOUBLE PRECISION FF (6) (input) : table of the function at six values of variable s stored 
in the array XSQ(6) (see below the commonblock M0NTAB1). 
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DOUBLE PRECISION VS(input) : input value of s. 

DOUBLE PRECISION FUN (output) : output value of the function. 

Calling by: MONDATA 



SUBROUTINE M0NLI2D (FDT , XI , X2 , Nl , N2 , XV , FUN) 



Purpose: Linear 2D interpolation from the table of any function. 

DOUBLE PRECISION FDT(N1 ,N2) (input) : Nl xN2 dimensional table of values for the in- 
put function. 

DOUBLE PRECISION XI (Nl) ,X2(N2) (input) : arrays for discrete points corresponding to 

the values of the input function. 
INTEGER Nl ,N2 (input) : dimensions of the 2D interpolation table. 
DOUBLE PRECISION XV (2) (input) : input values for two variables of the function. 
DOUBLE PRECISION FUN(output) : output value of the function. 
Calling by: MONDATA 



SUBROUTINE M0NLI4D (FDT , XI , X2 , X3 , X4 , Nl , N2 , N3 , N4 , XV , FUN) 



Purpose: Linear 4D interpolation from the table of any function. 

DOUBLE PRECISION FDT(N1 ,N2,N3,N4) (input) : NlxN2xN3xN4 dimensional table of 

values for the input function. 
DOUBLE PRECISION XI (Nl) ,X2(N2) ,X3(N3) ,X4(N4) (input) : arrays for discrete points 

corresponding to the values of the input function. 
INTEGER Nl ,N2 ,N3 ,N4 (input) : dimensions of the 4D interpolation table. 
DOUBLE PRECISION XV (4) (input) : input values for four variables of the function. 
DOUBLE PRECISION FUN (output) : output value of the function. 
Calling by: MONDATA 



SUBROUTINE M0NIN2D (FF , XI , X2 , Nl , N2 , FF1 , FF2 , FF3) 



Purpose: calculations of additional integrated tables used in the generation subroutine 

M0NGE2D,M0NG2D4. 
DOUBLE PRECISION FF(N1 ,N2) (input) : input table of 2D function. 
DOUBLE PRECISION XI (Nl) ,X2(N2) (input) : arrays for discrete points corresponding to 

the values of the input function. 
INTEGER Nl ,N2 (input) : dimensions of the 2D interpolation table. 

DOUBLE PRECISION FFl(Nl) ,FF2(N1) , FF3 (2 ,N1 ,N2) (output) : generated auxiliary ta- 
bles. 

Calling by: MONDATA, M0NIN4D 



SUBROUTINE M0NIN4D (FF , XI , X2 , Nl , N2 , FF1 , FF2 , FF3) 



Purpose: calculations of additional integrated tables used in the generation subroutine 
M0NG2D4. 

DOUBLE PRECISION FF(N1 ,N1 ,N2,N2) (input) : input table of 4D function. 

DOUBLE PRECISION XI (Nl) ,X2(N2) (input) : arrays for discrete points corresponding to 

the values of the input function. 
INTEGER Nl,N2(input) : dimensions of the 4D interpolation table. 

DOUBLE PRECISION FF1 (2 ,N1 , Nl ,N2) ,FF2 (2 ,N1 , Nl ,N2) , FF3 (8 , Nl , Nl , N2 , N2) (output) : 
generated auxiliary tables. 
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Calling by: MONDATA 



SUBROUTINE MONDATA 



Purpose: to read tables for absorptive corrections and for pRn form factors from the 
external files Spi_l, Sro_l, Sa2_l, S2pi_l, S2ro_l, S2a2_l, FFpi_l, FFro_l, 
FFa2_l. These tables are used for calculation of the differential cross sections 
for SCE and DCE reactions at given energy (defined by parameter M0NPAR(3)) 
by interpolation methods. 

Status of call: is called if M0NPAR(7)=1 or M0NPAR(8)=1. 

Calling by: MONINIT 

Calling of: MONCUBI, M0NLI2D, M0NLI4D, M0NIN2D, M0NIN4D 

3.3 Main Functions 



DOUBLE PRECISION FUNCTION MONCSEC(KP,KR) 



Purpose: to give the value of the total cross section of SCE (pp — > nX) or DCE (pp — > 
nXn) reaction for the given reggeon exchange at the c.m.s. energy defined 
by parameter M0NPAR(3) for the model defined by parameters M0NPAR(4) and 
M0NPAR(5). 

INTEGER KP {input) : single or double exchange 

= 1 : for SCE cross section 
= 2 : for DCE cross section 

INTEGER KR {input ) : type of the reggeon exchange 

= 1 : for SCE define 7r + exchange, for DCE n + ii + one. 
= 2 : for SCE p+ exchange, for DCE n+p + . 
= 3 : for SCE exchange, for DCE n+a^,. 

Calling by: MONINIT 



DOUBLE PRECISION FUNCTION MONCSCE (NO , XI , QT) 



Purpose: to give the value of cross section of SCE (pp — > nX) reaction for the given 
reggeon exchange at given £ n of neutron and Qt of reggeon at the c.m.s. en- 
ergy defined by parameter M0NPAR(3) for the model defined by parameters 
M0NPAR(4) and M0NPAR(5). 

INTEGER NO {input) : type of the reggeon exchange 

= 1 : for 7r + exchange. 
= 2 : for p + exchange. 
= 3 : for a\ exchange. 

DOUBLE PRECISION XI {input ) : £ n = |pbea '"~ p " 1 , relative momentum loss of the neu- 

Pbeam 

tron. 

DOUBLE PRECISION QT {input) : Q t , transverse momentum of the exchange reggeon. 
Calling by: MONCDCE, MONDATA 



DOUBLE PRECISION FUNCTION MONCDCE (NO , XII ,XI2 , QT1 , QT2) 
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Purpose: to give the value of cross section of the DCE (pp — > nXn) reaction for the given 
reggeon exchange at given of neutron and Q t ' of reggeons at the c.m.s. 
energy defined by parameter M0NPAR(3) for the model defined by parameters 
M0NPAR(4) and M0NPAR(5). 

INTEGER NO (input ) : type of the reggeon exchange 

= 1 : for 7r + 7r + exchange. 
= 2 : for 7r + p + exchange. 
= 3 : for 7T + a2 exchange. 



DOUBLE PRECISION XI (.input ) 
trons. 

DOUBLE PRECISION QT (input ) 
Calling by: MONCDCE, MONDATA 



Sri 



1,2 



"beam ^ n 

T71 



relative momentum loss of the neu- 



t>be 



Ql' 2 , transverse momentum of the exchange reggeons. 



DOUBLE PRECISION FUNCTION MONCSRP (NO, NCSMOD, SVAR) 



Purpose: to give the value of the total reggeon-proton cross section 
INTEGER NO (input ) : type of the reggeon exchange 

= 1 : for 7r + exchange. 
= 2 : for p + exchange. 
= 3 : for a~2 exchange. 

INTEGER NCSMOD (input ) : type of model for the reggeon-proton cross section calcula- 
tion 

= 1 : Donnachie-Landshoff parametrization [22]. 
= 2 : COMPETE parametrization [23]. 
= 3 : Bourrely-Soffer-Wu parametrization [24J. 
= 4 : Godizov-Petrov parametrization [25] . 

DOUBLE PRECISION SVAR (input ) : invariant mass of reggeon-proton system 
Calling by: MONCSRR, MONCSCE, MONCDCE 



DOUBLE PRECISION FUNCTION MONCSRR (NO , NCSMOD , SVAR) 



Purpose: to give the value of the total reggeon-reggeon cross section 
INTEGER NO (input ) : type of the reggeon-reggeon exchange 

= 1 : for 7r + 7r + exchange. 
= 2 : for ir + p + exchange. 
= 3 : for 7T + a^" exchange. 

INTEGER NCSMOD (input ) : type of model for the reggeon-reggeon cross section calcu- 
lation 

= 1 : Donnachie-Landshoff parametrization [22]. 
= 2 : COMPETE parametrization [23]. 
= 3 : Bourrely-Soffer-Wu parametrization [24J. 
= 4 : Godizov-Petrov parametrization [25] . 

DOUBLE PRECISION SVAR (input ) : invariant mass of reggeon-reggeon system 
Calling by: MONDATA, MONCDCE 
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3.4 Main Commonb locks and Parameters 



PARAMETER (MXGLPAR=200) 
REAL MONPAR 

COMMON/MONGLPA/ MONPAR (MXGLPAR) 



Purpose: to give access to the main MoNCHERswitches and parameters 



MONPAR(l) : number of events for the generation. 

MONPAR (2) : switch for LHE output. 

= : LHE output is switched off. 

= 1 : LHE output is switched on. 

MONPAR (3) : pp centre mass energy, in GeV, (from 900 to 14000 GeV). 

MONPAR (4) : kod of model for pR and RR interaction. 

= 1 : Donnachie-Landshoff parametrization [22]. 

= 2 : COMPETE parametrization [23]. 

= 3 : Bourrely-Soffer-Wu parametrization [24J. 

= 4 : Godizov-Petrov parametrization [25] . 

MONPAR (5) : kod of model for absorptive corrections. 

= 1 : 3 Pomerons eikonal model |13j . 

MONPAR (6) : type of exchange reggeon. 

= 1 : for SCE define 7r + exchange, for DCE 7r + 7r + one. 

= 2 : for SCE p+ exchange, for DCE tt+ p + . 

= 3 : for SCE exchange, for DCE % + 02- 

MONPAR (7) : switch for SCE generation. 

= : SCE is switched off. 

= 1 : SCE is switched on. 

M0NPAR(8) : switch for DCE generation. 

= : DCE is switched off. 

= 1 : DCE is switched on. 



Note 1: if M0NPAR(7)=0 and M0NPAR(8)=0, minimum bias events are generated 
by Pythia. 

Note 2: in the present version of Moncher, v. 1.1, the simultaneous generation 
of SCE and DCE is impossible. 



DOUBLE PRECISION S 

INTEGER NM0DPP , NM0DRR , ITYPR 

C0MM0N/M0NTAB0/S , NM0DPP , NM0DRR , ITYPR 



Purpose: to give access to some important Moncher parameters. 



S : pp c.m.s. energy, in GeV. 

NMODPP : kod of model for absorptive corrections. 

= 1 : 3 Pomerons eikonal model [l~3] . 

NMODRR : kod of model for pR and RR interaction. 

= 1 : Donnachie-Landshoff parametrization [22]. 

= 2 : COMPETE parametrization [23]. 

= 3 : Bourrely-Soffer-Wu parametrization [24j. 

= 4 : Godizov-Petrov parametrization [25] . 

ITYPR : type of exchange reggeon. 

= 1 : for SCE define 7r + exchange, for DCE 7r + 7r + one. 
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= 2 : for SCE p+ exchange, for DCE n+p+. 
= 3 : for SCE a 2 exchange, for DCE 7r + a^. 



DOUBLE PRECISION XSQ , PI , MPI , MP , MN , MRHO , MA2 
C0MM0N/M0NTAB1/XSQ (6), PI, MPI, MP, MN, MRHO , MA2 



Purpose: to give access to some important Moncher parameters. 

XSQ : six values of y/s for the interpolation subroutine MONCUBI. 

PI : 3.141592653589793D0 

MPI : pion mass. 

MP : proton mass. 

MN : neutron mass. 

MRHO : p meson mass. 

MA2 : a 2 meson mass. 



DOUBLE PRECISION XIMIN , XIMAX , QTMIN , QTMAX 
C0MM0N/M0NTAB2/XIMIN , XIMAX , QTMIN , QTMAX 



Purpose: to give access to some important Moncher parameters. 

XIMIN : minimal value of the variable £. 

XIMAX : maximal value of the variable £. 

QTMIN : minimal value of the variable |q| (transverse momentum of the neutron). 

QTMAX : maximal value of the variable |q| (transverse momentum of the neutron). 



DOUBLE PRECISION API , ARHO , AA2 , R2PI , R2RH0 , R2A2 
C0MM0N/M0NTAB3/ API , ARHO , AA2 , R2PI , R2RH0 , R2A2 



Purpose: to give access to some important Moncher parameters. 

API : slope of the pion regge trajectory. 

ARHO : slope of the p meson regge trajectory. 

AA2 : slope of the a 2 meson regge trajectory. 

R2PI : slope of the exponent in the residue of the pion trajectory. 

R2RH0 : slope of the exponent in the residue of the p meson trajectory. 

R2A2 : slope of the exponent in the residue of the a 2 meson trajectory. 



DOUBLE PRECISION GPI , GRHO , GA2 , SIGRSQ , KARHO , KAA2 
C0MM0N/M0NTAB4/GPI , GRHO , GA2 , SIGRSQ , KARHO , KAA2 



Purpose: to give access to some important Moncher parameters. 

GPI, GRHO, GA2 : constants Gl +pn /{8ir), G 2 p+pn /{&n) and G 2 a +J{%n). 

SIGRSQ : \r] R \ 2 . 
KARHO, KAA2 : k p , K a2 . 
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DOUBLE PRECISION SSPI , SSRHO , SSA2 , SDPIS , SDPIA , 
& SDRHOS , SDRHOA , SDA2S , SDA2A , FFDP I , FFDRHO , FFDA2 

C0MM0N/M0NDGET/SSPI(6,53,41) , SSRHO (6 , 53 , 51) , 
& SSA2(6,53, 51) , SDPIS (6, 10,8,9,8) ,SDPIA(6, 10,8,9,8) , 
& SDRHOS (6, 10,8,9,8) ,SDRH0A(6, 10,8,9,8) ,SDA2S (6, 10,8,9,8) , 
& SDA2A (6, 10,8,9,8) ,FFDPI(6,60, 16) , FFDRHO (6, 60, 16) , 
& FFDA2(6,60,16) 



Purpose: to give access to the input tables. 



DOUBLE PRECISION XSSPI ,XSSRHO ,XSSA2 ,XSDPIS ,XSDPIA , 
& XSDRHOS , XSDRHOA , XSDA2S , XSDA2A , XFFDPI , XFFDRHO , XFFDA2 

COMMON/MONDFIX/XSSPI (53,41) ,XSSRH0(53,51) , 
& XSSA2(53,51) ,XSDPIS(10,8,9,8) ,XSDPIA(10,8,9,8) , 
& XSDRHOS (10, 8, 9, 8) ,XSDRH0A(10,8,9,8) ,XSDA2S(10,8,9,8) , 
& XSDA2A(10,8,9,8) , XFFDPI (60, 16) , XFFDRHO (60 , 16) , 
& XFFDA2(60,16) 



Purpose: to give access to the additional tables obtained from the input files. 



DOUBLE PRECISION SPI , SRHO , SA2 , DPI , DRHO , DA2 , FDPI , FDRHO , FDA2 
C0MM0N/M0NDM0D/SPI(41,41) ,SRH0(41,41) ,SA2(41,41) , 

& DPI (17, 17, 17, 17) ,DRH0(17,17,17,17) ,DA2(17, 17, 17, 17) , 

& FDPI (17, 17) ,FDRH0(17,17) ,FDA2(17,17) 



Purpose: to give access to the tables for 2D and 4D generations. 



DOUBLE PRECISION SPIX1 , SPIX2 , SPIXQ , SR0X1 , SR0X2 , SROXQ , 
& SA2X1 , SA2X2 , SA2XQ , DPIX1 , DPIX2 , DPIXX , 
& DR0X1 , DR0X2 , DROXX , DA2X1 , DA2X2 , DA2XX 

C0MM0N/M0NDGE1/SPIX1 (41) , SPIX2 (41) , SPIXQ (2 ,41 ,41) , 
& SR0XK41) ,SR0X2(41) ,SR0XQ(2,41 ,41) , 
& SA2XK41) ,SA2X2(41) , SA2XQ (2 ,41 ,41) , 
& DPIXK17) ,DPIX2(17) ,DPIXX(2,17,17) , 
& DR0XK17) ,DR0X2(17) , DROXX (2, 17, 17) , 
& DA2XK17) ,DA2X2(17) ,DA2XX(2, 17, 17) 



Purpose: to give access to the auxiliary tables for 2D and 4D generations. 



DOUBLE PRECISION DDPI1 ,DDPI2 ,DDPI3 ,DDR01 ,DDR02 ,DDR03 , 
& DDA21 , DDA22 , DDA23 

C0MM0N/M0NDGE2/ 
& DDPI1(4,17,17,17) ,DDPI2(4,17,17,17) ,DDPI3(8, 17, 17, 17, 17) , 
& DDR01(4,17,17,17) ,DDR02(4, 17, 17, 17) ,DDR03(8, 17, 17, 17, 17) , 
& DDA21(4,17,17,17) ,DDA22(4, 17, 17, 17) ,DDA23(8, 17, 17, 17, 17) 



Purpose: to give access to the auxiliary tables for 2D and 4D generations. 
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DOUBLE PRECISION VXIR,VFIS,VFIA,VQTR, 
& VXIRF , VFIF , VXI , VQT , SVXI , SVQT , DVXI , DVQT 

COMMON/MONDVAR/VXIR(10) ,VFIS(8) ,VFIA(8) ,VQTR(9) , 
& VXIRF (60) ,VFIF(16) , VXI (53) ,VQT(41) , 
& SVXI (41) ,SVQT(41) ,DVXI(17) ,DVQT(17) 



Purpose: to give access to the arrays of variables for the input and auxiliary tables. 
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4 Program Installation 

Some materials related to the Moncher physics and generator is the one found on the 
web page 



http : //rioutine . web . cern. ch/rioutine 



in the section " Generators" . To get the code of the generator one should download the 
file 



http : //rioutine .web . cern. ch/rioutine/gencode/moncherl . 1 . tar .gz 



The program is written essentially entirely in standard Fortran 77, and should run on any 
platform with such a compiler. 

The following installation procedure is suggested for the Linux users, it was tested 
with CERN SLC5. 

$ gunzip moncherl . 1 .tar .gz 
$ tar -cvf moncherl . 1 .tar 
$ cd moncher/1 .1.0 
$ Is 

Now you can see some files: 

README contains brief description of the files in the current directory; 
moncher . f is the code of the generator; 

moncher . par defines switch keys and parameters for the simulation; 

Spi_l Sro_l Sa2_l contain data for the calculations of absorptive corrections for SCE; 

Sp2i_l S2ro_l S2a2_l contain data for absorptive corrections for DCE; 

FFpi_l FFro_l FFa2_l contain data for form-factors; 

mkmoncher is the executable file to compile and link moncher. f; 

rmoncher is the executable file to run moncher created by mkmoncher. 

$ ./mkmoncher 

compiles moncher. f by #77 compiler and link the generator with Pythia 6.420 [F5] and 
some CERNLIB libraries. Then, created executable moncher should be run by 

$ ./rmoncher 

Result of the simulation should be the PYTHlAstandard listing of one generated event of 
the SCE reaction pp — > nX at c.m.s. energy 7 TeV. The listing should be printed on the 
screen. If you have passed successfully all above, get start with the next step. 
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5 Getting Started with the Simple Example 

The Simple Example could look as following: 
PROGRAM MAIN 

IMPLICIT DOUBLE PRECISION (A-H, 0-Z) 
IMPLICIT INTEGER (I-N) 
c... global MONCHER parameters 
INTEGER MXGLPAR 
REAL MONPAR 

PARAMETER (MXGLPAR=200) 
COMMON/MONGLPA/ MONPAR (MXGLPAR) 



c. . .initialization 

CALL M0NGIVE('M0NPAR(1)=1000') 
CALL MONGIVE ('MONPAR (2) =1') 
CALL MONGIVEC MONPAR (3) =7000') 
CALL MONGIVE ('MONPAR (4) =1') 
CALL MONGIVE ('MONPAR (5) =1') 
CALL MONGIVE ('MONPAR (6) =1') 
CALL MONGIVE ('MONPAR (7) =1') 
CALL MONGIVE ('MONPAR (8) =0') 
CALL MONGIVE ('MSEL=2') 
CALL MONINIT 

NT0T=M0NPAR ( 1 ) 
KLHE=M0NP AR ( 2 ) 



number of events 

switch for LHE saving 

pp centre mass energy in GeV 

code of model for pR/RR interaction 

code of model for absorption 

type of Reggeon 

switch for SCE generation 

switch for DCE generation 

pythia: mb+sd+dd+elastic+lowpt 



c . . .generation 

DO NEV=1,NT0T 
CALL MONEVEN 

IF(NEV.EQ.l) CALL PYLIST(l) 
CALL ANALYZER (I OUT) 

IF(KLHE.EQ.l.AND.IOUT.EQ.l) CALL MONUPEV 
ENDDO 

c... final statistics 
CALL PYSTAT(l) 

c... produce final Les Houches Event File. 
IF(KLHE.EQ.l) CALL PYLHEF 

STOP 
END 

First, we set some values for elements of array MONPAR which control a process of gen- 
eration. Then, we should initialize the generator calling MONINIT. In this example we are 
going to generate 1000 events of Single Pion Exchange, pp — > n{ii + p) — > nX, at c.m.s. en- 
ergy 7 TeV. The (vr + p) interaction is controlled by Pythia and it includes minimum bias, 
single and double diffraction, elastic scattering and low-pt scattering. Filling of MONPAR 
elements can be done also from the external file moncher.par. Subroutine M0NPARA call- 
ing by MONINIT checks the presence of the moncher.par in the current directory and, if 
it exists, reads parameters MONPAR, see chapter [6j 
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On the next step, we generate some number of events, defined by M0NPAR(2). Every 
event is generated by MONEVEN. User's subroutine ANALYZER (I OUT) is called after every 
event generation, analyses the event and sets some value to the integer variable IOUT If 
I OUT is equal to unity, we save this event in the LHE format using the subroutine MONUPEV. 

Here you can see example of the Simple Analyzer: 

SUBROUTINE ANALYZER (IOUT) 
IMPLICIT DOUBLE PRECISION (A-H, 0-Z) 
IMPLICIT INTEGER(I-N) 
C...HEPEVT commonblock. 

PARAMETER (NMXHEP=4000) 

COMMON/HEPEVT/NEVHEP , NHEP , ISTHEP (NMXHEP) , IDHEP (NMXHEP) , 
& JMOHEP (2 , NMXHEP) , JDAHEP (2 , NMXHEP) , PHEP (5 , NMXHEP) , VHEP (4 , NMXHEP) 
DOUBLE PRECISION PHEP, VHEP 
SAVE /HEPEVT/ 

c 

IOUT =0 
I SIGN =1 
NEUTR0NS=0 

c 

CALL PYHEPC(l) 

c 

DO 1=1, NHEP 
KP =IDHEP(I) 
ETA =PYP(I,19) 

IF (KP.EQ. 21 12. AND. DABS (ETA) .GE.8.5) THEN 
NEUTR0NS=NEUTR0NS+1 
ISIGN=ISIGN*ETA 
END IF 
ENDDO 

c 

IF (NEUTRONS . EQ . 2 . AND . ISIGN . LT . 0) I0UT=1 

c 

RETURN 
END 

In this example, we analyse all particles in the generated event and look for the neutrons 
(code 2112) in the region of pseudorapidity \rj\ > 8.5 (assumed acceptance of the neutron 
detector). If number of such neutrons is equal to 2 and they move in opposite directions, 
IOUT is set to unity. 

Finally, we print the PYTHIAstatistics by PYSTAT and produce the final LHE file which 
has the name moncher.lhe by default. 

This example has a concrete physical meaning. We have selected SCE events with 2 
leading neutrons moving in the opposite directions which imitate a DCE process. So, we 
have saved background for the DCE from the SCE. 
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6 Program Control Parameters 



All parameters that control the generation can be defined in the external file moncher . par. 
For example, the set of parameters for the generation of the S7rE process, described in 
the chapter [5l can look as follows: 



c MONCHER v. 1.1.0 card file 

c 

c MONCHER control keys 

c 

M0NPAR(1)=1000 ! number of events to generate 
c 

M0NPAR(2)=1 ! key for Les Houches data(l-save,0-no) 
c 

M0NPAR(3)=7000 ! pp centre mass energy in GeV (900 -> 14000) 
c 

M0NPAR(4)=1 ! code of model for pR and RR interaction 

c NM0DRR=1 -> Donnachie-Landshoff model (default) 

c NM0DRR=2 -> COMPETE (PDG) model 

c NM0DRR=3 -> Bourreli-Sopf er-Wu model 

c NM0DRR=4 -> Godizov-Petrov model 

c 

M0NPAR(5)=1 ! code of model for absorption 

c NM0DPP=1 3 IP eikonal model (default) 

c not now NM0DPP=2 -> Godizov-Petrov model 

c not now NM0DPP>2 -> other models . . . 

c 

M0NPAR(6)=1 ! type of Reggeon (l-pi+, 2-rho+, 3-a2+) 

c (for DCE only pi-pi, pi-rho and pi-a2 survive) 

c 

M0NPAR(7)=1 ! key for SCE generation 
c 

M0NPAR(8)=0 ! key for DCE generation 
c 

c PYTHIA control keys 

c 

cMSEL =0 ! full user control 

cMSUB(ll)=l ! f + f -> f + f (QCD) 

cMSUB(12)=l ! f + fbar -> f + fbar' 

cMSUB(13)=l ! f + fbar -> g + g 

cMSUB(28)=l ! f + g -> f + g 

cMSUB(53)=l ! g + g -> f + fbar 

cMSUB(68)=l ! g + g -> g + g 

cMSUB(91)=l ! Elastic scattering 

cMSUB(92)=l ! Single diffractive (AX) 

cMSUB(93)=l ! Single diffractive (XB) 

cMSUB(94)=l ! Double diffractive 

cMSUB(95)=l ! Low-pT scattering 

cMSEL =1 ! mb 

MSEL =2 ! mb+sd+dd+elastic+lowpt 
c 

MRPY(1)=12031967 ! start point of random number generator 
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Subroutine MONPARA reads lines from moncher.par. All lines begining with a letter 
"c" are ignored by the program, all others lines are processed by subroutine MONGIVE, 
which can recognize any variables from the MONCHERcommon block /MONGPGL/ and 
the PYTHIAcommon blocks /PYJETS/, /PYDAT1/ , /PYDAT2/, /PYDAT3/ , /PYDAT4/ , 
/PYDATR/ , /PYSUBS/, /PYPARS/, /PYINT1/, /PYINT2/, /PYINT3/ , /PYINT4/, 
/PYINT5/, /PYINT6/ , /PYINT7/, /PYINT8/, /PYMSSM/, /PYMSRV/, /PYTCSM/, 
/PYPUED/, (see [15]) . Parameters MONPARA are described in detail in the section |3l 
page IT51 

Using parameters from the common blocks listed above, one can define wide spectrum 
of SCE (M0NPAR(7)=1) and DCE ((M0NPAR(8) = l) processes or any processes existing in 
PYTHIA (if (M0NPAR(7)=0 and M0NPAR(8)=0). Some examples are described in the next 
chapter. 

7 Examples of the Moncher Processes. 



N 


Process 


Type of n + p 
interactions 


Picture of the process 


The 
Moncher 
parameters 


1 


pp — > nX 


minimum bias: 

7T + P — > X 


n * 

P , 


*^ P 
X 


M0NPAR(7)=1 
M0NPAR(8)=0 
MSEL=1 


2 


pp — > nn + p 


elastic scattering: 

7l + P — > 7l + p 


P . 

n+ 

P 


r^^^^ n 

p 


M0NPAR(7)=1 
M0NPAR(8)=0 
MSEL=0 
MSUB(91)=1 


3 


pp — > nXY 


double 
diffraction: 
n+p -> X + Y 


P^^ 
Y^^ 


^-^^^ n 
^P 


M0NPAR(7)=1 
M0NPAR(8)=0 
MSEL=0 
MSUB(94)=1 


4 


pp — >■ nXp 


single diffraction 
(7r + dissociation): 

7l + P — > X + p 


P^^ 

P--^~ 


^-^^ n 

p 


M0NPAR(7)=1 
M0NPAR(8)=0 
MSEL=0 
MSUB(92)=1 


5 


pp — > nXii + 


single diffraction 
(p dissociation): 
n + p -> X + 7T+ 


P 

X 


71+ 


M0NPAR(7)=1 
M0NPAR(8)=0 
MSEL=0 
MSUB(93)=1 



Table 4: Some S7rE processes which can be generated with Moncher. 



It was mentioned already in Chapter [3] that the Moncher generates pEn vertices 
and, then, IRp (for SIRE) or 1R1R (for DIRE) interactions are generated by Pythia. The 
type of these interactions can be controled by the PYTHlAparameters. We can define 
elastic or inelastic interactions, diffractive or non-diffractive processes, different types of 
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diffraction, hard scattering, etc. Some of the basic processes for S7rE and D7rE, which can 
be generated by the Moncher, are presented in the tables EH and [5] respectively. 

Let us consider one more simple example, how to generate process number 2 from 
Table HI This is a Single Pion Exchange with elastic scattering of the virtual pion by the 
proton of the beam. This reaction, pp — > mr + p, has very clear signature: neutron, proton, 
single 7r + meson and nothing else in the final state. Initial particles are scattered at very 
small angles and, thereof, there are no any detector signals in the region of pseudorapidity 

< 7. An experimental possibility of such measurements has been analysed in Ref. [8] 
with prereleased version of Moncher. 

File moncher . par with parameters for the generation of pp — > nir + p can look as follows: 



M0NPAR(1)=1 

M0NPAR(2)=0 

M0NPAR(3)=7000 

M0NPAR(4)=1 

M0NPAR(5)=1 

M0NPAR(6)=1 

M0NPAR(7)=1 

M0NPAR(8)=0 

MSEL =0 

MSUB(91)=1 



number of events to generate 

key for Les Houches data(l-save ,0-no) 

pp centre mass energy in GeV (900 -> 14000) 

code of model for pR and RR interaction 

code of model for absorption 

type of Reggeon (l-pi+, 2-rho+, 3-a2+) 

key for SCE generation 

key for DCE generation 

full user control 

elastic scattering 



Parameter M0NPAR(7)=1 defines the generation of the SIRE process. Exchange reggeon is a 
pion (M0NPAR(6)=1). PYTHIAparameters MSEL=0 and MSUB(91)=1 set elastic 7c + p scatter- 
ing. Parameter M0NPAR(4)=1 sets Donnachie-Landshoff parametrization for 7i + p interac- 
tion, see subsection 12. 1.21 Parameter M0NPAR(5)=1 specifies 3 Pomeron model for absorp- 
tive correcttions, see subsection 12. 1.1 [ Parameters M0NPAR(1)=1 and M0NPAR(3)=7000 set 
the generation of 1 event at 7 TeV pp c.m.s. energy. We don't ask to save any events 
(M0NPAR(2)=0) and the only result of the generation is the PYTHlAlisting of the generated 
event: 

Event listing (summary) 
I particle/ jet KS KF orig p_x p_y p_z E 



m 



1 !p+! 21 2212 0.000 0.000 3500.000 3500.000 0.938 

2 !p+! 21 2212 0.000 0.000-3500.000 3500.000 0.938 



3 


nO 


1 


2112 


2 





114 





216- 


-2296 


804 


2296 


804 





940 


4 


!pi+! 


21 


211 


2 


-0 


114 


-0 


216- 


-1203 


196 


1203 


196 





140 


5 


!p+! 


21 


2212 


3 


-0 


019 


-0 


001 


3500 


000 


3500 


000 





938 


6 


!pi+! 


21 


211 


4 


-0 


095 


-0 


215- 


-1203 


196 


1203 


196 





140 


7 


P+ 


1 


2212 


5 


-0 


019 


-0 


001 


3500 


000 


3500 


000 





938 


8 


pi+ 


1 


211 


6 


-0 


095 


-0 


215- 


-1203 


196 


1203 


196 





140 






sum: 


2.00 







000 





000 





000 


7000 


000 


7000 


000 



In this listing lines 1 and 2 correspond to the protons of the beams. Lines 3, 7 and 8 
relate to the neutron, proton and pion, respectively, in the final state of the reaction. 
The proton is deflected at angle ^5.5xl0 -6 rad., neutron and pion are scattered in the 
direction opposite to proton, as it is shown on the diagram of the process in the table HI 
with polar angles ^10 -4 and ^2xl0 -4 rad. 
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N 


Process 


Type of 7T + 7r + 
interactions 


Picture of the process 


The 

MONCHER 

parameters 


1 


pp — > nXn 


minimum bias: 

7T + 7T+ ->■ X 


P - 

»+y 


^ n 
* n 


M0NPAR(7)=0 
M0NPAR(8)=1 
MSEL=1 


2 


pp — > nn + 7i + n 


elastic scattering: 

7T + 7T + — > 7T + 7T + 


P , 

;r+ n: 
n 


^-^^^ n 

J 7T+ 
P 


M0NPAR(7)=0 
M0NPAR(8)=1 
MSEL=0 
MSUB(91)=1 


3 


pp — > nXYn 


double 
diffraction: 

7T + 7T+ -> X + Y 


P „ 

n 


n 

7t+ * 

P 


M0NPAR(7)=0 

1 1U111 XXX b \ 1 / \j 

M0NPAR(8)=1 

MSEL=0 

MSUB(94)=1 


4 


pp — > nX7r + n 


single diffraction: 

7T + 7T+ X + 7T + 


P 

n 


n 

.7T+ • 
P 


M0NPAR(7)=0 
M0NPAR(8)=1 
MSEL=0 
MSUB(92)=1 

or 

MSUB(93)=1 



Table 5: Some D7rE processes which can be generated with MONCHER. 
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Index of Subprograms and Common Block Variables 
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subroutine 
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subroutine 


ma 
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subroutine 


nn 


MONMBDF 


subroutine 


QH 


MONEVEN 


subroutine 


QH 


MONSPEG 


subroutine 


QH 


MONDPEG 


subroutine 


QJ 


MONSPEM 


subroutine 


d2 


MONDPEM 


subroutine 


d2 


MONSHPY 


subroutine 


QJ 


MONGIVE 


subroutine 


QJ 


MONUPEV 


subroutine 


US] 


MONUPIN 


subroutine 


m 


M0NGE2D 


subroutine 


M 


M0NG2D4 


subroutine 


M 


MONCUBI 


subroutine 


M 


M0NLI2D 


subroutine 


m 


M0NLI4D 


subroutine 


QJ 


M0NIN2D 


subroutine 


m 


M0NIN4D 


subroutine 


m 


MONDATA 


subroutine 


M 


MONCSEC 


function 




M 


MONCSCE 


function 




M 


MONCDCE 


function 




m 


MONCSRP 


function 






MONCSRR 


function 






MONGLPA 


common 


block 


M 


MONTABO 


common 


block 
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M0NTAB1 


common 


block 
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M0NTAB2 


common 


block 
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M0NTAB3 


common 


block 


m 


M0NTAB4 


common 


block 


m 


MONDGET 


common 


block 


M 


MONDFIX 


common 


block 


m 


MONDMOD 


common 


block 


ED] 


M0NDGE1 


common 


block 


ED] 


M0NDGE2 


common 


block 


Eoj 


MONDVAR 


common 


block 


m 


MONPAR in /MONGLPA/ 


CD 


S in /MONTABO/ 




da 


NMODPP in /MONTABO/ 


QJ 


NMODRR in /MONTABO/ 


us 


ITYPR in /MONTABO/ 
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